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DENSITY GRADIENTS AND INTERNAL DUST IN THE ORION NEBULA 



El resumen sera traducido al espanol por los editores. The ionization structure of the Orion nebula 
can be described as a skin-like ionization structure on the surface of a dense cloud. We propose that a steep 
density stratification, increasing as a powerlaw (n oc a; -2 ) function of distance x from the ionization front, 
exhibits properties which agree with our long-slit spectrum of the Orion nebula. For instance, there exist a 
unicity relation between both the H/3 surface brightness or the ionization front [Sli] density, and the scale 
L of the powerlaw, where L is the distance between the ionization front and the onset of the density near 
the exciting star. Internal dust is required to obtain a simultaneous acceptable fit of both the [S n] density 
and then H/3 suface brightness observations. Nebular models containing small dust grains provide a better fit 
than large grains. The line ratio gradients observed along the slit are qualitatively reproduced by our density 
stratified models assuming a stellar temperature of 38000 K. Collisional deexcitation appears to be responsible 
for half of the gradient observed in the [Nii]A5755/[Nii]A6583 temperature sensitive ratio. We propose that 
the empirical relationship found by Wen & O'dell between density and stellar distance may possibly be caused 
by a power-law density stratification. 



The ionization structure of the Orion nebula can be described as a skin-like ionization structure on the surface 
of a dense cloud. We propose that a steep density stratification, increasing as a powerlaw (n oc x~ 2 ) function 
of distance x from the ionization front, exhibits properties which agree with our long-slit spectrum of the Orion 
nebula. For instance, there exist a unicity relation between both the H/3 surface brightness or the ionization 
front [S il] density, and the scale L of the powerlaw, where L is the distance between the ionization front and 
the onset of the density near the exciting star. Internal dust is required to obtain a simultaneous acceptable 
fit of both the [Sn] density and then H/3 suface brightness observations. Nebular models containing small dust 
grains provide a better fit than large grains. The line ratio gradients observed along the slit are qualitatively 
reproduced by our density stratified models assuming a stellar temperature of 38000 K. Collisional deexcitation 
appears to be responsible for half of the gradient observed in the [N n]A5755/[N n]A6583 temperature sensitive 
ratio. We propose that the empirical relationship found by Wen & O'dell between density and stellar distance 
may possibly be caused by a power-law density stratification. 
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RESUMEN 



ABSTRACT 



It was recognized early on by Munch (1958) and 
Wurm (1961) that the Orion nebula is better de- 
scribed as a photoionized skin on the surface of a 
dense cloud than by a uniformly filled gas struc- 
ture. Furthermore, kinematics of the ionized gas (see 
review by O'Dell 1994) indicates the gas is flowing 
away from the dense ionizing front on the surface of 
the molecular cloud. Although earlier studies pro- 
posed an exponential function to describe the den- 
sity behaviour from such a flow (Tenorio-Tagle 1979; 
Yorke 1986), more recent theoretical and observa- 
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1. INTRODUCTION 



tional works (Franco, Tenorio-Tagle & Bodenheimer 
1989, 1990; Franco et al 2000a,b) indicate that ex- 
panding H II regions are well described by a power 
law density stratification n oc x~@ with exponents 
larger than (3 = 1.5. This is qualitatively consistent 
with the results obtained from hydrodynamic calcu- 
lations of photoevaporative flows off proplyd disks by 
Richling & Yorke (2000) or off spherical surfaces by 
Bertoldi (1989) and Bertoldi & Draine (1996). The 
current paper looks at how the observed gradient of 
the [S n] density in Orion can be accounted for using 
a simple powerlaw density model and how well the 
surface brightness compares with the long-slit obser- 
vations presented here. The aim is to explore the 
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benefits of a powerlaw density stratification in the 
context of a determination of the shape of the ion- 
ization front (hereafter IF) using a surface brightness 
map in Ha, a task first performed by Wen & O'Dcll 
(1995) for the two cases of isochoric and exponential 
density behaviour. 



2. OBSERVATIONS AND DATA REDUCTION 

We have taken long-slit spectra of the cen- 
tral region of Orion at the Observatorio Astroffsico 
Guillermo Haro, Cananea, on the 21st of Decem- 
ber 1998 using a 1024 x 1024 TEK CCD detector 
mounted on a Boiler & Chivens spectrograph. We 
used a grating of 300 1/mm and covered the wave- 
length range 3600-5200 A, 5500-7100 A, and 8000- 
9600 A in three settings. The corresponding expo- 
sure times were 90, 60 and 300 seconds, respectively. 
A slit width of 2.5" was used. The effective spectral 
resolution was « 6 A. The slit was aligned West-East 
and off-centered by 30" to the South of 8 1 C Ori as 
illustrated in Fig. [I]. Arc spectra were taken as well 
as standard star exposures. Data reduction (bias 
subtraction, flat fielding, wavelength and flux cal- 
ibration) was performed with IRAF. Due to non- 
photometric conditions and partial cloud coverage, 
we could not derive a reliable absolute flux calibra- 
tion. We hence used the calibration of Rodriguez 
(1999), who had also observed with the slit aligned 
East- West with her regions A-l and A-3 common to 
ours. We tied our H/3 fluxes to the mean surface 
brightness within her region A-l. In any event, the 
adopted Rodriguez surface brightness value is only 
25% higher than our initial calibration. Line fluxes 
were reddening corrected using the Balmer decre- 
ment and the galactic reddening law (Seaton 1979). 
Further details on the data can be found in Gonzalez- 
Gomez (1999). 

3. THE PHOTOIONIZATION MODELS 

We adopt a similar approach to that of Bincttc & 
Raga (1990, hereafter BR90) who studied the prop- 
erties of powerlaw densities in photoionized slabs (see 
also Williams 1992). However, we adapted their sim- 
ple slab geometry to the more appropriate spherical 
geometry for reasons given below in Sect. 3.2. An- 



other paper which considers powerlaw density strat- 
ification in photoionization models is that of Sankrit 
& Hester (2000) (see also Hester et al. 1996). 

The emission line measures in Orion are domi- 
nated by the gas near the IF as a result of the strong 
density gradient towards it. In the work of Wen 
k O'Dell (1995, hereafter W095), the authors ap- 
proximated the front structure as gas segments (or 



Fig. 1. Image illustrating the position of our long-slit 
on an R-band image of Orion. The position of the excit- 
ing star 8 1 C Ori (whose image is saturated) is indicated. 
The R-band filter encompasses Ha and the red [N n] lines 
which account for most of the nebular emission visible on 
the image. 



columns) aligned towards our line of sight with a 
projected area on the sky of 2" x 2". In their fi- 
nal best model, the densities within the segments 
declined exponentially towards the observer rather 
than radially towards 6 1 C Ori. Using the observed 
distribution of Ha surface brightness and [Sll] gas 
density, W095 were able to solve for the IF distance 
(along our line of sight) for the whole nebula. Such a 
distribution of distances is termed the "shape" of the 
IF. The gas segments of their model did not extend 
up to the exciting star, a characteristic also shared 
by the isobaric slab model of Baldwin et al. (1991, 
hereafter BF91). In such a scheme, the segments 
are independent and a density must be specified for 
each segment. In this Paper, we will explore density 
distributions which are exclusively function of radius 
from 6 1 C Ori without assuming any significant gas 
cavity near 8 1 C Ori. 

3.1. The radial density gradient 

We adopt a similar notation to BR90 in which 
the density profile (with an origin lying inside the 
dense cloud) is described as follows: 



n(x) = n ;x < L, 

n(x) =0 x > L , 



(1) 
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Fig. 2. Diagram illustrating qualitatively our powerlaw 
density stratification as a function of x. The origin of the 
spatial x axis lies at a depth L from the nebular boundary 
and the density at that position (at x = L) is labelled 
n (see Equation Q). Note that L will be allowed to vary 
with angular direction from 8 1 C Ori. The top axis is the 
radial distance r from the star 8 1 C Ori, with r increasing 
towards the left. The IF (or inner Stromgren boundary) 
is situated at X s (or equivalently at R s ). The electronic 
density inferred from the [S n] doublet corresponds to a 
density near the IF. The hollow gas cavity surrounding 
the star has a size r . In our models, the adopted value 
for r is much smaller than depicted above. 

where n(x) is the gas density, j3 the powerlaw index 
(n oc x~ l3 ) 1 n a the boundary density of the nebula 
towards the exciting star, and x the variable repre- 
senting the distance from the inner density discon- 
tinuity at x = 0. Note the reversal of axis relative 
to traditional studies: the ionizing radiation is en- 
tering at the cloud boundary x — L and is absorbed 
inward up to the inner Stromgren boundary at X s 
(see Fig. ||). The steep density gradient near the ori- 
gin ensures that the photons are absorbed before the 
density discontinuity, hence X s > 0. 

In the case of Orion, the spherical geometry is 
necessary in order to take into account the geomet- 
rical dilution of the ionizing radiation across the ra- 
dial density structure. One can define the variable, 
r, as the distance from the central ionizing star; that 
is r = L — x + r a where r Q is the radius of the 
central hollow cavity and n the gas density at r Q . 
When the ionizing photon luminosity Qh is suffi- 
ciently high (or the density n low), the Stromgren 
radius, R s — L — X s + r Q , is near the origin of the 
x axis, at a distance ~ L + r a from the exciting star 
(since X s is negligible). This is referred by BR90 as 
the "strong gradient" regime. In the converse case, 
for very low ionizing photon luminosities Qh (or for 
high boundary densities n ), R s shrinks [X s becomes 



non-negligible) and becomes of order r Q + small frac- 
tion of L, a condition which BR90 labeled "weak 
gradient" because the nebula's thickness is smaller 
than the density gradient's scale L. All the quanti- 
ties introduced above are illustrated in Fig. [|. 

The "strong gradient" regime is the appropri- 
ate one for Orion since the nebula is traditionally 
approximated by a skin-like structure on the face 
of a dense molecular cloud (BR90; O'Dell 1994). 
One advantage of the powerlaw stratification when 
the regime of "strong gradient" prevails is that we 
know before hand the approximate radius of the IF 
(< L + r a ). In the case of an exponential strati- 
fication, which is much shallower than a powerlaw, 
it requires solving for the ionizing radiation transfer 
inside the nebula in order to predict the IF depth. 

For definiteness, we will consider the case of a 
powerlaw index j3 = 2. Very similar nebular line 
ratios of [O n] , [O in] , etc. were obtained with indices 
of j9 = 1.5 and 1.0, provided we adjusted n D and L 
in such a way that both the ionized depth R s and 
the [S n] density remained the same as that obtained 
with the (3 = 2 case. We consider that our results do 
not depend on the particular choice of j3 made here. 

We found no justification for having an arbitrary 
large region devoid of gas around 9 1 C Ori and, for 
definiteness, have set r Q = 0.007pc < 0.05R S , a size 
so small that its precise value has no impact on the 
calculations. 

3.2. The geometry 

The geometry of the IF is highly structured and 
complex (c.f. W095) and cannot be modeled di- 
rectly by a code like MAPPINGS ic which is designed 
to handle simple nebular geometries. Our approach 
has been to compute a sequence of spherical mod- 
els of different scale L (c.f. Equation |l|) in order 
to simulate small areas of the IF located at differ- 
ent distances from the exciting star. This piece- 
meal approach (a different model for each position 
along the slit) allowed us to reconstruct the IF con- 
tribution through each aperture element. To account 
for the projection of each slit aperture element onto 
the calculated Stromgren spheres, we carried out 
in MAPPINGS ic the integral of line emissivities over 
spatial limits corresponding to the appropriate aper- 
ture projection onto the calculated spheres. 

Our aim is to model only the background 
IF (the component B in BF91) which is much 
brighter than the foreground IF (the component A 
in BF91). We therefore halved the fluxes computed 
by MAPPINGS Ic to reduce the geometry to that of 
an hemisphere rather than a sphere. Furthermore, 
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to remove the small contribution due to component 
A, we multiplied the observed H/3 surface brightness 
by |, the same value as in BF91. 

To sum up, the different IF positions covered by 
the slit are associated to a sequence of models which 
differ only by their value of the scale L. Since our slit 
was offset to the South by 30" from 1 C Ori (extend- 
ing mostly eastward), it did not sample gas compo- 
nents which overlaid along common radial lines from 
the exciting star. 

3.3. Metallicity, dust and Qh 

Dust mixed with nebular g CIS, 3b possibility stud- 
ied by BF91, can play an important role in the 
structure of the Orion nebula. We define /i£> as 
the amount of dust (by mass) present in the gas 
in units of the solar neighborhood value. The code 
MAPPINGS ic includes the effects of internal dust as 
described in Binette et al. (1993). We will consider 
two extinction curves which present quite a different 
behavior in the far-UV, both calculated by Martin & 
Rouleau (1991): a) the orion curv eQ used by BF91 
is characterized by a population of relatively large 
grain sizes covering the range 0.03-0.25 /im, b) the 
noiro curve, first introduced by Magris, Binette & 
Martin (1993), is characterized by a population of 
small grain sizes covering the range 0.005-0.03 /im. 
A possible justification is that photoerosion, which is 
significant in photoionized gas, could have the effect 
of reducing the mean grain size to that characterizing 
the noiro curve. Interestingly, a distribution of small 
grains may also play an important role in resolving 
the temperature problem in Hn regions (Stasihska 
& Szczerba 2001). Note that the solar neighborhood 
extinction curve as calculated by Martin & Rouleau 
(1991) is composed of grain sizes covering the full 
range 0.005-0.25 /im (i.e. ism = orion + noiro). The 
advantage of using the noiro and the orion curves is 
that they provide a simple characterization of fairly 
opposite cases regarding extinction properties in the 
UV. 

We adopt the same gas phase abundances as 
those determined by BF91 (mean metallicity at 
0.7ZQ) and the stellar atmosphere models of Hum- 
mer & Mihalas (1970) with T eff = 38000 K, using 
the scheme proposed by Shields & Searle (1978) to 
interpolate in effective temperature and metallicity 
between atmosphere models. The choice of temper- 
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ature is discussed in Sect. 4.3 



As argued by W095, the observed recombination 
lines surface brightness includes a reflected back- 

1 The orion extinction curve is certainly indicated for the 
cold gas in front of the Orion nebula (inside the "Lid" ) but is 
not necessarily appropriate for the warm ionized gas. 
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Fig. 3. Behavior of the density inferred from the [S n] 
lines as a function of the scale L employed in different 
photoionization calculations (each model is represented 
by a filled square). The long-dashed and short-dashed 
line represent two dust-free sequences of models in which 
the boundary density is n a = 150 and 2000 cm -3 , respec- 
tively. The dotted line and the solid line represent dusty 
models (all with a boundary density n = 150) in which 
the orion (/id = 0.8) or the noiro (fin = 0.4) extinction 
curve was employed, respectively. 



component due to scattering from the background 
dusty molecular cloud. To correct for this, as in 
W095 we multiplied the observed H/3 surface bright- 
nesses by | and adopted a reduced stellar photon 



luminosity of Qh = 10 49 photons s 
(see discussion in W095). 



in the models 



4. RESULTS OF THE MODELS 

Our main proposition is that only one basic vari- 
able of the nebular geometry, the scale L, varies be- 
hind the slit aperture. We will show that the obser- 
vations we have taken, albeit very limited, are not 
inconsistent with this simple picture. 

4.1. Properties of a powerlaw stratification 

For a given dust content fip and dust extinction 
curve, we proceed to show that for a given value of 
n , there exists a unicity relation between the [Sn] 
density (which is representative of the IF electronic 
density) and the scale of the gradient L. This is 
shown in Fig. || where we plot the [S Ii] density as a 
function of the powerlaw scale L employed in differ- 
ent photoionization calculations. 
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Fig. 4. Figure illustrating the sky projection of the West- 
East aligned slit onto the nebular IF. The distance be- 
tween 1 C Ori and the IF is R B ■ The substellar distance 
is d (see Sect. |4.2| ) . There is a further offset of the slit to 
the South of flTj Ori of 0.073 pc (not shown). 

In the dust-free case, the long-dashed line cor- 
responds to the sequence of models with n = 
150 cm -3 while the short-dashed line corresponds 
to models with a higher boundary density n a — 
2000cm -3 . For the denser model, the density gra- 
dient is much too shallow since we get a R s of only 
0.25 pc in the model with L = 0.7 pc. Clearly, the pa- 
rameter n a plays an important role. Its value can be 
deduced by requiring that the best model sequence 
fits simultaneously the observed profiles along the 
slit in H/3 brightness and in [S n] density. 

In Fig. ||, the dotted line represents the case of 
using the orion extinction curve with /i£>=0.8 (this 
sequence of models have the same dust character- 
istics as those assumed by BF91). The use of the 
smaller grains noiro curve (solid line) on the other 
hand causes a considerable change in the inferred 
position of the IF even though these models con- 
tain only half the mass of grains present in the orion 
curve models (dotted line). The reduction of the IF 
radius is due to the amount of ionizing photons ab- 
sorbed by the dust grains. 

In their determination of the IF shape, W095 
found a very interesting albeit unexplained relation- 
ship between the distance to 9 1 C Ori and the gas 
density (their Fig. 6). They also found that the 
closer the IF is to 9 1 C Ori, the higher the Ha bright - 
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Fig. 5. Behavior of the [S n] density after projection onto 
our long- slit (D axis). The slit was aligned EW. D — 
corresponds to the pixel nearest to 8 1 C Ori with D in- 
creasing towards the East, crossing the Bright Bar near 
D = 0.25. The long and short dash line models represent 
the dust-free case with n a = 150 and 2000 cm -3 , respec- 
tively. The dotted and the solid lines represent models 
with internal dust using the orion and noiro curves, re- 
spectively (both with n — 150). The two dot-dash lines 
represent the noiro dusty model using an boundary den- 
sity of n a — 50 (short-dash dotted line) and 500 cm -3 
(long-dash dotted line), respectively. 

ness appears to be. As shown in Fig. |^, a powerlaw 
stratification establishes a relationship between the 
scale L and the [S n] density near the IF. A similar 
relationship holds between the H/3 surface bright- 
ness and the [S II] density. Our proposed powerlaw 
stratification is therefore apt in providing a credi- 
ble explanation to the intriguing [S n] density-radius 
correlation discovered by W095. 

4.2. Comparison with our long-slit data 

Comparison of the models with our long-slit data 
will be done considering both the [S n] density and 
the H/3 surface brightness. One difficulty resides in 
relating the local IF radius, R s , with position, D, 
along our slit. We have proceeded as follows. We 
have not resolved the inverse problem of finding the 
distance of the IF at each observed position along D, 
Our main goal consists rather in exploring whether 
a powerlaw density is generally consistent with the 
data at hand. With this perspective in mind, we 
settled for the following simple relation between D 
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Fig. 6. Behavior of the H/3 surface brightness in powerlaw 
density stratified models after projection onto our long- 
slit. The long-dash and short-dash lines represent dust- 
free models with n — 150 and 2000 cm -3 , respectively. 
The dotted and solid line (both with n a = 150 cm -3 ) 
are models with internal dust using the orion and noiro 
curves, respectively. The two dot-dash lines represent the 
noiro dusty model with boundary densities of n = 50 
(short-dash dotted line) and 500 cm -3 (long-dash dotted 
line), respectively. The observed H/3 surface brightnesses 
plotted here have been reduced to corr ect for the contri- 
bution of the scattered light (Sect. 3.3) and for the con- 
tribution of the foreground component (Sect. 3.2), as in 
W095 and BF91, respectively. 



and R s : D — yjR 2 s - 0.073 2 pc - d where 0.073 pc 
is the slit offset of 30" to the South, R s the true IF 
distance to 9 1 C Ori for a given slit position D in par- 
sec projected on the sky. D = corresponds to the 
pixel nearest to 9 1 C Ori with D increasing towards 
the East. The additional parameter d (in pc) was 
obtained by iteration, using as constraint that our 
best model should fit as closely as possible both the 
observed [Sil] densities and the H/3 surface bright- 
nesses. If the IF lied on a perfect plane perpendicular 
to our line-of-sight, d would correspond to the sub- 
stellar distance (O'Dell 1994) and it would have to 
be subtracted in quadrature from r (as for the slit 
offset). This resulted, however, in an unsatisfactory 
fit for all the models tried. A plausible explanation 
can be found in the 3D-map of the IF shape of W095 
which shows that the IF plane is tilted away from us 
at the corresponding position of our slit. The above 
equation, in which d is subtracted linearly, qualita- 
tively describes a tilted plane. A cartoon describing 
the proposed IF geometry is shown in Fig. [|. 

In Fig. |, we compare our stratified models with 
the [Sn] densities observed along the slit. For sim- 
plicity, all models cover the same range in L, from 
0.22 to 0.66 pc, which is the appropriate range for 
the favored solid line model. This has the advantage 
of clarifying the effect of varying only one parameter 
at a time (either n Q or ud)- The surface bright- 
nesses along the slit are shown in Fig. ^ and can be 
compared with the various model sequences. The 
two Figures were drawn using d = 0.26 which is the 
value required by the favored model (solid line). 

We emphasize that no dust-free model could pro- 
vide an acceptable fit to the data in both figures || 
and ||, whatever value was used for d. In the case of 
dusty models using the alternative extinction curve 
orion, the fit was never satisfactory in both Figures 
simultaneously. In effect, if we increase d to 0.35, 
this shifts all model sequences to the left by 0.09 pc 
in both Figures, thereby resulting in a good fit of the 
orion model in Fig |B| but at the expense of an un- 
satisfactory fit in Fig ^. Increasing the dust content 
brings the dusty models down but it is not possible 
in the case of the orion since ud = 0.8 is already 
the maximum allowed given the subsolar metallicity 
of Orion. The preeminence of small grains, on the 
other hand, presents the advantage that less dust is 
necessary (ud = 0.4 for the noiro curve). This sce- 
nario makes sense in the event that dust grains have 
been progressively photoeroded. 

In conclusion, even when allowing d to vary 
freely, the Orion model consisting of boundary den- 
sity n a = 150 cm -3 and containing small grains 
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{noiro with /Xd= 0.4) provides the optimized set of 
physical conditions^ favored by our model for the 
following reasons: 

1 . dust-free models are unsuccessful in reproducing 
the observed behavior of the H/3 surface bright- 
ness of Fig. |6| (adopting a larger value for Qh 
would worsen the discrepancy), 

2. for a given dust content and size distribution, 
the determination of the boundary density n 
is straightforward since the curves using differ- 
ent densities lie relatively far apart in the [S n] 
density diagram of Fig. || (for instance, compare 
the position of the two dash-dotted lines corre- 
sponding to sequences with boundary densities 
of n — 50 or 500 cm -3 ), 

3. models favoring the preeminence of small grains 
(e.g. the noiro extinction) provide a simultane- 
ous fit to both data sets presented in Figs || and 
||, respectively. 

The applicability of the proposed model does not 
include the so-called "Bright Bar" (> 0.25 pc) be- 
yond which the H/3 brightness is slightly rising. Ac- 
counting for this rise in our simple picture would 
require a gas structure which lies closer to 9 1 C Ori 
than the calculated IF model above. Alternatively, 
the rise may be due to some sort of limb brightening. 
Dopita et al. (1974) have proposed for instance that 
the Bright Bar correspond to the IF seen edge-on. 
This is qualitatively corroborated by the IF distance 
map produced by W095, which shows a sharp ridge 
at the position of the Bar. 

4.3. Gas temperature and line ratios 

Although the current work does not attempt to 
provide a fit of all the observed line ratios, it is 
worth checking whether the proposed powerlaw den- 
sity stratification results or not in ratios similar to 
the observed. Because the ionization parameter in 
our models is not a free parameter, line ratios be- 
come an interesting consistency check. Furthermore, 
it is worth testing whether the models can reproduce 
the observed line ratio trends with slit position. 

In a similar fashion to BF91, one important 
parameter of our models is the stellar tempera- 
ture which directly affects the excitation line ratio 
[O in] /H/3. We found that reproducing satisfactorily 
this ratio requires a stellar temperature of 38000 K 

2 Interestingly, the apparent reddening manifested in the 
emergent line ratios due to internal dust is very small. In 
effect, the emergent Balmer decrement for the small-grain and 
the dust-free models are 2.97 and 2.88 respectively 



when using the stellar atmosphere models]] of Hum- 
mer & Mihalas (1970). This temperature is consis- 
tent with an evolved star of spectral type 07 (c.f. 
W095). Our line ratio calculations of [Oinl/H/3 is 
represented by the solid line in panel d of Fig. 7] using 
the same parameters as before, that is a boundary 
density n — 150 and an extinction curve consisting 
of small grains (noiro) with /irj = 0.4. We discuss 
below possible explanations for the solid line lying 
slightly below the data. In any event, it is encourag- 
ing that the models show the same gradient as the 
data. 

As for the temperature indicators 
Rnii = [Nii]A5755/[Nii]A6583 and R om = 
[Olll]A4363/[Oin]A5007, a gradient is also present 
in the data as shown by panels b and c (BF91 
reported also a gradient of R NII although for a dif- 
ferent slit alignment). According to the models, the 
gradient in the ratio Rom is due to a temperature 
gradient (AT ~ 500 K) while in the case of Nitrogen 
it is due to both a density gradient (via collisional 
deexcitation, c.f. BF91) and a temperature gradient. 
Although the stratified models can successfully fit 
the observed trends in panels b and c, the predicted 
ratios lie somewhat lower than observed, more so 
for Rom- Lowering the metallicity would increase 
Rom but at the expense of having the [O in] /H/3 
ratio becoming larger than observed. 

It is well known that photoionization models do 
not explain why the temperatures derived from col- 
lisionally excited lines differ from that of recombi- 
nation lines (see review by Peimbert 1995). In the 
case of Orion, Esteban et al. (1998) found that the 
abundance of the +2 ion derived from the O n re- 
combination lines are incompatible with the same 
ion abundance derived from the collisionally excited 
[O m] A5007 when a common temperature is adopted. 
Deriving a unique ionic abundance was found to re- 
quire different mean temperatures. These can be 
justified if there exist high amplitude temperature 
fluctuations in the nebula (Peimbert 1967) which far 
exceed the amplitude arising from traditional pho- 
toionization models (Kingdon & Ferland 1995, 1998; 



3 More recent stellar atmosphere models exist (Hubeny & 
Lanz 1995; Schaerer & de Koter 1997; Hillier & Miller 1998; 
Pauldrach et al. 2001) but they differ considerably according 
to whether they include or not the effect of stellar winds [c.f. 
Fig. 5 in Sankrit & Hester (2000)]. Until the issue of wind 
loss is resolved, all stellar temperatures inferred from pho- 
toionization models will be of limited accuracy. The models of 
Hummer & Mihalas were used above for their convenience as 
they can be interpolated to any arbitrary temperature within 
MAPPINGS Ic. Our results do not depend on which family of 
models is used provided the stellar temperature is adjusted to 
reproduce the line ratios well. 
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Fig. 7. The dots represent the observed line ratios. Panel a: the line ratio [0 n] AA3727/[0 ffl]A5007 which is a good 
indicator of the degree of ionization, panel b: the Nitrogen temperature sensitive ratio [N n]A5755/[N ll]A6583, panel c: 
the Oxygen temperature sensitive ratio [O ill] A4363/[0 iii]A5007, panel d: the [O ill] A5007/H/? A4861 ratio. The scale 
of the 3 panels b, c and d is the same and extend over 0.42 dex. The models (see Sect. 4.3): the solid line shows the 
behavior of line ratios along the slit for our most successful model (n — 150 and the noiro curve with fio = 0.4). The 
long-dashed line was calculated assuming temperature fluctuations with t 2 = 0.01. The dotted line was calculated using 
half the diffuse field from the "outward-only" approximation. The spike in the data near D — —0.05 pc in panels b and 
c is believed to be due to a foreground proplyd. 
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Perez 1997). The mean square fluctuation ampli- 
tudes that Esteban et al. inferred, lie in the range 
t 2 = 0.02-0.028. The effect of temperature fluctu- 
ations on line cmissivitics were incorporated in the 
code MAPPINGS ic along the scheme described in Bi- 
nette & Luridiana (2000). The long-dashed lines 
in Fig. [7] represent such a model calculated with 
t 2 = 0.01. As expected, Rnii and, more so, Rom 
are both significantly increased and can now repro- 
duce satisfactorily the data of panels b and c. The 
calculated [Om]/H/3 ratio in panel a, however, lies 
now too high and the stellar temperature would have 
to be reduced in order to produce a self-consistent 
model in all panels. Given our current poor knowl- 
edge about the nature and cause of the postulated 
temperature fluctuations, we have not carried out 
further the iterative process. 

The diffuse ionization radiation field plays a non- 
negligible impact on the gas temperature. In our 
models, we use the "outward-only" approximation 
to treat the transfer of the diffuse radiation. This 
is equivalent to assuming a single direction for the 
propagation of the diffuse field, a technique based 
on the approximation that the diffuse field is soft 
and does not therefore travel far from the point of 
emission. To test whether the integration along the 
direction of increasing density may cause an overes- 
timation of the intensity of the diffuse radiation, we 
have calculated models in which we reduced by half 
the strength of the diffuse field (adopting t 2 = 0). 
Such models are represented by the dotted-line in 
Fig. 0. We can see a significant shift which is in the 
same direction as that caused by temperature fluc- 
tuations with t 2 ~ 0.01. The dotted line gives us 
a useful upper limit on the line ratio uncertainties 
resulting from a possible inadequacy of the outward- 
only approximation. 

Panel a is a plot of the [Oii]/[Om] ratio. The 
models can reproduce the observed points only 
within a restricted slit range. No obvious explana- 
tion can be found for the discrepancy outside it. It 
possibly reflects the limitations of our simplified ge- 
ometry as discussed in Sect. 3.2 or the larger errors 
affecting the UV [0 11] lines. Overall, the fit of the 
gradients using our simple models appears satisfac- 
tory in all the panels of Fig. ^. 

5. CONCLUSIONS 

The method used by W095 to determine the IF 
shape requires, at each point of the nebula, knowl- 
edge of both the Ha surface brightness and the IF 
density (using the [S11] or [Oil] doublet). For a 
powerlaw density stratification as proposed in this 



Paper, the density across the nebula is not a free pa- 
rameter since it is uniquely determined by the local 
powerlaw scale L. This opens the possibility that 
using only the Ha (or H/3) surface brightness map, 
one could in principle determine the IF shape of the 
whole nebula. For this reason, further studies are 
warranted in order to: a) extend the comparison of 
powerlaw density models to a more extensive set of 
observations of Orion and b) work out a scheme to 
solve for the inverse problem (as was done by W095 
in the exponential case) in order to determine the IF 
shape assuming the proposed powerlaw density gra- 
dient and taking into account the effects of internal 
dust. 
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